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ABSTRACT 


A comparison is made between the relative levels of aircraft interior 
noise related to structureborne and airborne paths for the same propeller 
source. A simple, but physically meaningful, model of the structure treats 
the fuselage interior as a rectangular cavity with five rigid walls. The 
sixth wall, the fuselage sidewall, is a stiffened panel. The wing is modelled 
as a simple beam carried into the fuselage by a large discrete stiffener repre- 
senting the carry-through structure. The fuselage interior is represented by 
analytically-derived acoustic cavity modes and the entire structure is repre- 
sented by structural modes derived from a finite element model. The noise 
source for structureborne noise is the unsteady lift generation on the wing 
due to the rotating trailing vortex system of the propeller. The airborne 
noise source is the acoustic field created by a propeller model consistent 
with the vortex representation. Comparisons are made on the basis of interior 
noise over a range of propeller rotational frequencies at a fixed thrust. The 
measure of noise level is based on nonresonant response to eliminate the large 
variations associated with unquant if table damping levels. It is found that 
the relative importance of the structural and airborne paths is highly depen- 
dent on the structural parameters. Over the range of parameters considered in 
this study it is found that the structureborne contribution can vary from very 
insignificant to nearly equivalent to the airborne levels. 



INTRODUCTION 


The question addressed here is the relative importance of structureborne 
and airborne paths for propeller noise in aircraft interiors. Of particular 
interest are wing-mounted engines driving multiple -blade propfans. 

A significant design consideration in the installation of advanced propfan 
propulsion systems on transport category aircraft is the interior noise level 
within the passenger compartment. The extrapolation of experience gained with 
conventional propeller -driven aircraft and the current generation of turboprop 
aircraft indicates that in the case of the advanced, highly loaded, multiple - 
blade propellers currently proposed the interior noise levels may exceed com- 
fort levels and may even exceed hearing damage levels. It is therefore impor- 
tant to understand not only the noise generating mechanism, but also the paths 
by which the noise is transmitted into the aircraft interior. 

Advanced turbofan engines are likely to be tractor or pusher wing-mounted 
installations or pusher aft mounted installations. The tractor wing-mounted 
engine, typical of current turboprop aircraft, has the most significant struc- 
ture-borne noise implications because of the interaction of the propeller 
trailing vortex system with the wing structure. It is this configuration 
which is considered here. 

For several years there has been considerable speculation on the relative 
importance of airborne and structureborne paths for the transmission of pro- 
peller noise into the fuselage. The interaction of the propeller rotating 
trailing vortex system with the wing structure creates an oscillatory forcing 
function on the wing which can be transmitted into the fuselage interior 
through the wing and fuselage structural coupling to the cavity. The airborne 
path describes the mechanism by which the propeller radiated noise field is 
transmitted through the fuselage sidewall to the fuselage cavity. 

No significant data base exists on the relative importance of the struc- 
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tureborne and airborne paths. Measurements to isolate the two sources on an 
unmodified aircraft are not possible with current experimental methods. Iden- 
tification of the two sources will probably require physical isolation, for 
example, breaking the structural path by disconnecting the wing from the fuse 
lage or breaking the airborne path by placing a barrier between the propeller 
and the fuselage. Tests of this type are major undertakings and have not to 
date been accomplished to the point of providing useable data. 

Very little work has been done on the modelling of the s true tureborne 
noise paths. In contrast, fairly extensive models have been developed for 
airborne noise [for example; 1,2,3], Metcalf and Mayes [4] addressed the 
question of the possibility of a significant contribution of structureborne 
noise to the overall interior noise levels based on tests which broke the air 
borne noise path by wrapping the fuselage of the test aircraft with an acous- 
tical barrier material. In a related study, Unruh [5] examined the structure 
borne noise path which transmits engine vibrations in a single engine general 
aviation aircraft by an experimental technique in which the engine was iso- 
lated from the fuselage. 

Until recently no substantial attempt has been made to produce an analyti 
cal model of the structureborne path for propeller noise. Junger, Garrelick, 
Martinez and Cole [6] appear to be the first to consider the introduction of 
propeller wake disturbances into the wing structure and subsequent wave propa 
gation into the fuselage structure. Their model was based on a Green's func- 
tion approach treating the wing structure as an acoustic wave guide, as dis- 
cussed by Junger and Feit [7] . 

In the present study a model is created which can be used to compare 
structureborne and airborne noise levels in the fuselage cavity. The mecha- 
nism for the introduction of propeller noise into the fuselage is basically 
the same as that used by Junger, Garrelick, Martinez and Cole [6]. The air- 
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borne noise source level is also modelled* here using a finite element propel- 
ler radiation model which predicts the radiated acoustic field of the propel- 
ler [8]. This model can be used in the near field of the propeller. For rel- 
atively small fuselage -propeller tip clearances the near field is important. 

The structural-acoustic system is based on a simple geometry which approx- 
imates the essential features of the wing- fuselage structure and the cabin 
interior. The interior is modelled as a rectangular cavity. The sidewall is 
modelled as a flat stiffened panel. The wing is modelled as a beam which is 
carried into the fuselage structure by a heavy discrete stiffener. The wing- 
fuselage is represented by a finite element procedure from which the free 
vibration frequencies and normal modes are determined. The cavity is repre- 
sented in terms of its hardwall acoustic natural frequencies and normal modes. 
For the rectangular cavity, or for the logical extension to a cylindrical 
cavity, the acoustic modes can be represented analytically. A coupling 
procedure discussed by Dowell, Gorman, and Smith [9] is used to construct the 
complete system model. 

The interior noise levels can be calculated in response to the structure- 
borne and airborne noise sources. The approach used here is to consider a 
particular propeller geometry operating at a fixed thrust and a given forward 
speed. This fully defines the airborne noise levels and the strength of the 
propeller vortex system for the s true tureborne source. The sound pressure 
level at a reference point in the cavity is calculated for a range of propel- 
ler rotational speeds. For a fixed thrust, this means that blade twist and 
loading are changed. The comparison of interior noise levels is based on the 
level of the nonresonant response to eliminate uncertainties created by 
unquant if iable levels of structural damping. 

Heitman and Mixson [10] have found that interior acoustic absorption may 
play an important role in controlling the levels of airborne interior noise. 


3 


In the present formulation a model of this type of acoustic treatment has been 
included. The results of computations including interior absorption will be 
reported later in this report. 

THE MODEL FOR THE STRUCTUREBORNE NOISE SOURCE 

The mechanism for the generation of structureborne noise is the interac- 
tion of the propeller tip vortex system with the wing. This interaction 
creates a spanwise variation of lift which varies periodically in time. 

Figure 1 shows, in simplified form, the trailing vortex system from a two- 
bladed propeller. The vortex system, created by the lift on the propeller 
blades, is swept around with the rotating propeller. The velocity field of 
each vortex modifies the downwash on the wing and, therefore, locally creates 
a periodic lift variation. 

The velocity field behind the propeller is very complicated, but the 
important features of the interaction with the wing can be modelled by assum- 
ing that there is associated with each propeller a vortex, the circulation of 
which is determined by the lift per unit span at the propeller tip. The bound 
vorticity on the propeller and the vortex sheet immediately behind the pro- 
peller are neglected so that only the fully developed rolled up tip vortex 
system is considered to be important. This description of the vortex system 
is consistent with the model of Junger, Garrelick, Martinez, and Cole [6]. 

The physical picture is that of an isolated vortex trailing from each 
propeller blade and rotating with the propeller. In an axis system centered 
on the vortex the induced velocity field is 

V = (r/ 27 rh) e^- 

where F is the circulation strength of the vortex, h is the distance from the 
vortex center to the point at which the velocity V is calculated, and e^- is 
a unit vector normal to the line between the vortex center and the point at 
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which V is calculated (i.e. a unit vector tangent to the circular streamline 
of the vortex) . 

With this model it is possible to represent the normal wash on the plane 
Z - 0 (in which the axis of the propeller lies) in four regions on the wing: 

0 < f} < 1: The region behind the propeller disk and outboard of the propeller 

hub 


"--S 


cosiNnt 


rf > 1: Outboard of the propeller disk 


w 




^-(iN+l) cosiNnt 


( 1 ) 


( 2 ) 


-1 < f? < 0: The region behind the propeller disk and inboard of the 

propeller hub 


w = cosiNflt 

i=l 


ry < -1: Inboard of the propeller disk 


where 




(.l)iN+l |^|-(iN-l) cosiNnt 


r? = (Y-Yo)/R 

and 


Yq = Spanwise location of the propeller axis 
Y = Spanwise location of the normal wash point 
R = Propeller radius 


(3) 


(4) 
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Q - Propeller rotational speed ( radians /sec. ) 

N - Number of propeller blades 

w — Normal wash (positive in the direction of positive z, upward) 

and 

w » (T/2nR) w (5) 

Figure 2 shows the spanwise variation of normal wash amplitude. 

The tip vortex strength T is determined from the lift per unit span at the 

tip 

it - pvr (6) 

where 

it lift per unit span at the tip 
V = helical velocity at the tip 
p = air density 

The circulation can thus be written as 

r = + (nR/u)2] V2) ( 7 ) 

where U is the forward velocity. The parameter it depends on the propeller 
model used. Junger, Garrelick, Martinez, and Cole [6] use a triangular load- 
ing which when combined with a given total lift, number of blades, forward 
flight dynamic pressure, propeller radius, and rotational speed, will define 
the entire propeller loading, and in particular, the tip loading. This type 
of loading is consistent with the actual loading on conventional propellers 
[ 11 ]. 

The resulting lift per unit span is 

i = (1/2) pu2(l + J-2)c^c^ at-(r/R) (8) 

and the resulting thrust loading and torque loading per unit span are 

t = (1/2) U^(1 + J a^(r/R)cos6 (9) 

m = 1/2 pV^a + J-^)C^Cj, af.(r/R)sinfi (10) 
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where the angle the angle between the helical velocity and the propeller 
disk, is defined by 

cosy3 - (r/JR)[l + (r/R)2j*2] -1/2 (H) 

sirv9 = [1 + (r/R)2j-2]-l/2 (12) 

The parameter J is defined by 

J » U/QR (13) 

and is proportional to the advance ratio used in the propeller literature. 
Also defined here are 

= propeller tip chord 
= tip angle of attack 

Cq =» propeller lift curve slope at the time 
The value of the tip angle of attack is determined by the total thrust where 
T = |[(1+ j2)l/2 . j2i,,[j-l + (1 + J*2)l/2]| (14) 

where 

it “ (1/2) pV^a + J'hc^Co «t (15) 

Equations (7) and (14) are sufficient to determine F, given the propeller 
thrust, geometry, and rotational speed. Equations (9) and (10) can then be 
used to compute the propeller loading. The loading is required in the calcu- 
lation of the airborne sound pressure levels. 

The normal wash distribution given by Eqs. (l)-(5) locally creates an 
effective angle of attack of the wing given by 

= W/U (16) 

which creates a lift force on the wing. In this analysis we use the quasi- 
static assumption and calculate the lift according to a strip theory repre- 
sentation as 

= (1/2) ^u2c„Co (17) 

where 
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Cn - 2 jt, the theoretical two dimensional lift curve slope for incompres- 
sible, steady flow 
- lift per unit 
« wing chord 

Junger, Garrelick, Martinez, and Cole [6] use a correction for unsteady 
flow and compressibility. However, this is not considered to be essential in 
the present analysis since such a refinement is certainly lost in other 
approximations which are required. Since the lift distribution is directly 
proportional to the downwash amplitude distribution, it follows that Figure 2 
also can be viewed as a distribution of the oscillatory lift distribution on 
the wing. 

An examination of Eqs . (l)-(5) and Eqs . (16) and (17) shows that for a 
specific propeller rotational speed fl, the structureborne source frequencies 
are harmonics of NQ, where N is the number of propeller blades. These are 
exactly the harmonics of the airborne noise source. In making computations of 
interior noise related to the structureborne source it has been found that the 
interior noise levels are well defined with two or three harmonics being com- 
puted. 


THE MODEL FOR AIRBORNE NOISE 

In order to compare the relative importance of airborne and structural 
paths for interior noise, it is necessary to have a model for airborne noise 
which is consistent with the model for the structureborne noise source dis- 
cussed above. Such a model has been created by Eversman and Steck [8] . In 
this model the propeller is represented by a thrust loading and torque loading 
distribution in the propeller disk. The loading which is used is given by 
Eqs. (9) and (10) and is therefore consistent with the structureborne noise 
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source. The thrust and torque loading become the volumetric force distribu- 
tion in a finite element formulation of the convected wave equation. The 
finite element solution is carried out in cylindrical coordinates for each 
harmonic of the blade passage frequency. A wave envelope scheme is used to 
represent the radiation boundary condition both with and without the forward 
flight effect[12]. 

If t(r) and m(r) are the thrust and torque components of the propeller 
blade loading per unit span, given by Eqs. (9) and (10), then the propeller 
can be replaced by an equivalent volumetric force distribution in the propel- 
ler plane given by [8] 

00 

fx “ «(x) E (18) 


E 


where 


_ .(pc2)-l[t(r)/a](Nn/27r) [sin(iNnt)/iNfi] (20) 

= (pc2)-l[m(r)/a](Nn/2jr)[sin(2Nnt)/iNQ] (21) 

and ^NQ. The parameter a is the projection of the blade chord on the 

propeller disk given by 

a = c costp (22) 

where c is the local blade chord, <p is the blade twist measuring the angle of 
the propeller chord plane with respect to the propeller disk plane, and r is 
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time required for the propeller to sweep by a point on the propeller disk, 

r - a(r)/Qr (23) 

Equations (18) and (19) show that the volumetric force distribution is 
periodic with harmonics of the fundamental frequency NQ. Each temporal har- 
monic is associated with angular harmonics of the fundamental angular wave 
length 27 t/N. It is important to note that the source harmonic frequencies for 
the fixed rotational speed Q are the same as for the structureborne source. 

Figure 3 shows a typical propeller noise radiation pattern in the form of 
contours of constant acoustic pressure magnitude. The contours shown are in a 
plane of constant ^ in a cylindrical coordinate system with the symmetry axis 
along the propeller axis and with the thrust direction to the right. The con- 
tours shown here are for a six-bladed propeller with no forward flight effect. 
Figure 4 is a cartesian plot of the directivity of the same propeller on a 
line two propeller radii from the propeller axis. This is the location which 
is chosen to represent the noise levels to which the cabin sidewall is sub- 
jected. The curve of Figure 4 is normalized to 100 dB, but the absolute level 
is also computed so that the maximum sound pressure level on the reference 
line is known. This maximum pressure, which can easily be in the range of 
105-135 dB, is used as the airborne noise source. In the present study no 
effort is made to model the actual distribution of sound pressure level on the 
cabin sidewall and it is assumed that the sound pressure level is constant 
over the sidewall and is normally incident. 


STRUCTURAL MODEL AND COUPLING TO THE CAVITY 
The structural model attempts to represent the essential features of the 
structureborne path of a wing-fuselage system. Both the wing and the fuselage 
are modelled with finite elements. This is schematically depicted in Figure 
5. For simplicity, the fuselage is represented by a rectangular cavity with 
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all sides rigid except for the side adjacent to the wing. The wing is repre- 
sented by a beam and the fuselage carry- through structure is represented by a 
heavy vertical stiffener, as shown. The fuselage sidewall is assumed to be 
clamped at its edges and the model is assumed to be symmetric with respect to 
the vertical plane formed by the wing and stiffener. 

Both the wing and carry- through stiffener are modelled by one -dimensional 
Hermitian beam elements in bending. Linear extensional properties are also 
included. Each node has three degrees of freedom: A transverse displacement, 

a slope, and an inplane axial displacement. The fuselage sidewall is modelled 
by a series of four-noded sixteen-degree-of- freedom Hermitian flat plate ele- 
ments [13]. The variables at each node are (w, w^^, Wy, w^y) where w is the 
transverse displacement and subscripts denote the partial derivative of w with 
respect to that variable. Figure 6 shows the discretization of the portion of 
the fuselage panel that was modelled. Symmetry about the vertical stiffener 
is assumed so that the stiffener is shown at the right boundary of the fuse- 
lage panel. In this analysis the fuselage panel is assiamed to be stiffened by 
stringers and frames. However, because these stiffeners are not as large as 
the structural connection for the wing to the fuselage, they are modelled as 
smeared stiffeners [14], In Figure 6, the spacing of the vertical frames and 

/s. 

horizontal stringers is d and i, respectively. The finite element model can 

r- 

easily be extended to include discrete frames and stiffeners since the ele- 
ments used throughout have enough degrees of freedom to enforce all the compa- 
tibility conditions. Thus, no growth in dimensionality would be required. In 
addition, the Bogner, Fox, Schmit [13] elements have been coded to be applica- 
ble to cylindrical shell representations for the sidewall. There are obvious 
possibilities for growth of the present model. 

The structural model has been generated with a program written specifi- 
cally for this problem in order to take full advantage of the coupling tech- 
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niques introduced. In this section an outline of the procedure is given. 

The wing is modelled using Hermitian bending and linear axial deflection. 
The nodal degrees of freedom are bending deflection, bending slope, and axial 
deflection. If there are wing elements, then the wing has 3N^ degrees of 
freedom (the degrees of freedom at the attachment point of the wing to the 
carry- through structure are assigned to the sidewall). These degrees of free- 
dom are designated (W^) . The fuselage sidewall (plate) has ^(Np^+ l)(Np^+ 1) 
degrees of freedom, when Np^ and Np^ are the number of elements lengthwise and 
widthwise, respectively. Each node has four degrees of freedom: w, w^, w^, 
and sidewall modal displacement vector is {Wp}. Note that all of 

the sidewall degrees of freedom are not active due to the boundary constraint. 
Note also that stringers, frames, and the carry- through structure will not 
introduce new degrees of freedom. 

The dynamic equations of motion for the wing and sidewall are written as 


■ 

" 
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Qp 
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(24) 


The force vector {Q^} is generated from the propeller source model and {Qp} is 
generated from the cavity acoustic pressure acting on the sidewall. The vec- 
tor {Q^} cannot be partitioned in terms of wing deflections only, since the 
wing attachement point degrees of freedom are also assigned to the sidewall. 
The vector {Qp} is partitioned in terms of the sidewall deflections only. A 
free vibration analysis, based on Eq. (24), is formulated as 
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A subspace interation eigenvalue routine was used to solve the eigenvalue 
problem. This yields a sequence of structural eigenvalues and a sequence 

of structural eigenvectors {Wg)j_, which, when truncated in a suitable way, can 
be used to generate a modal matrix 


[^1 = 




(Ws)i {Ws) 2 . 



(26) 


This modal matrix is partitioned into wing and plate contributions . An 
eigenvector expansion of Eq. (24) yields 


^ 1 r 



Msii (q) + Cg 

_ ^ _ 

ii (q) + 

Ksii 

L ^ 


(q) - ['h/ » [^pl^{Qp) (27) 


where 



t^]T [M] [^] 


(28) 


is the diagonal generalized mass matrix for the structure and [M] is the ori- 
ginal finite element mass matrix. In the calculation of eigenvalues, [M] is 
normalized so that - 1, and An equivalent viscous modal 

damping Cg^j - 2f£U)£ is included at this stage. The coupling from the cavity 
to the plate in the structural equations is through the new generalized force 
vector 

(Qsc) “ t^p]'^{Qp) 

The form of this vector can be obtained by returning to the definition of 
(Qp). The virtual work of the cavity pressure on the sidewall is 


5W - 



s 


p(x,z)6Wp(x,z)dS 
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where p(x,z) is the cavity acoustic pressure, 5wp(x,y) is the virtual change 
in the transverse deflection of the plate, and 3Qg is the sidewall surface 
area. In terms of the implicit global shape matrix for the finite element 
representation of the plate 


6wp(x,z) - [Np(x,z)] {5Wp) 


so that 


= {5Wp) 


/ 


[Np]^p(x,z)dS 




The cavity acoustic pressure is given in terras of the cavity modes 


p(x,z) = [^]{a} 


where the modal matrix [(f)] is a suitably truncated sequence of the cavity mode 
shapes evaluated on the sidewall and (a) is a vector of modal amplitude coef- 
ficients . 

With these observations, the generalized force vector which couples the 
cavity to the structure is 


{Q 


sc 


} = t^p 



{a} = [D]{a) 


s 

The integration is broken down into subdomains which are the finite elements 
on the plate surface and the shape function matrix [N] is explicitly defined 
within these subdomains by the 16 degree -of- freedom plate element shape func- 
tions. The coupling matrix is denoted as [D] and has Ng rows and columns, 
where Ng is the number of retained structural generalized coordinates and N^ 
is the number of retained cavity modes . 
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CAVITY MODEL AND COUPLTNG TO THE PLATE 
For the present study the cavity Is rectangular in geometry with dimen- 
sions a (length), b (height), and c (depth). The cavity is to be represented 
in terms of analytic acoustic modes calculated for a completely rigid cavity. 
These modes and the associated natural frequencies can be calculated directly 
from the three dimensional wave equation. The cavity frequencies are given by 


, m , n 


:o»r (£/ay 


(m/b)- 


(n/c)- 


where Cq is the speed of sound in the cavity. The frequencies are arranged in 
a sequence co^, i - 1, 2, ... N^, in ascending order. The eigenfunctions are 


^ i,m,n (x,y,z) ■= cos(ijrx/a)cos(mjry/b)cos(nffz/c) (30) 


where the cavity coordinate system is shown in Figure 5. The eigenfunctions 
are also ordered in a sequence ^^(x,y,z), i - 1, 2, . . . N^, corresponding to 
the sequence of eigenvalues. It is important to note that the lowest cavity 
frequency is wq q q - - 0 and the corresponding cavity mode, or eigen- 

function, is ^]^(x,y,z) “ 1, that is, uniform pressure. 

The cavity modes satisfy the orthogonality condition 

/ <}>i,<f>\ “ Mcii^iJ 

where 


and 



= abc , i - m = n - 0 
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^i>0, m-0, n-0 
Mcii - abc/2 , i i-0, m>0, n-0 
^i«0, m-0, n>0 

^i>0, m>0, n«0 
-= abc/4 i — 0, m>0, n>0 

L^>0, m«0, n>0 

= abc/8 , i>0, m>0, n>0 

In this study, is referred to as "acoustic generalized masses" for the 

cavity, although the units are those of volume. 

Coupling of the structure to the cavity is based on a weighted residuals 
procedure which reduces to the method discussed by Dowell, Gorman, and Smith 
[9]. 

The acoustic pressure within the cavity is governed by the wave equation 

v2p - (l/Co2) in n (31) 

and the boundary condition 

Vp • n “ -PqW^^. on dQ 

where Q is the volume to which Eq . (26) applies, and dQ is the surface of H, 
with outward normal n, w^^- is the outward particle acceleration at the cavity 
wall, and is the density of the air in the cavity. Boundary conditions of 
relevance include 
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( 32 ) 


V p • ^ — 0 on 
A>o(A/Po‘^)Pt 

where 3flj| represents hard walls, dOg the flexible sidewall, and 50 l 
represents acoustically lined area where the relation between acoustic pres- 
sure and particle velocity is 

Pt “ Z wt-t (33) 

and the nondimens ional normal incidence admittance is given by 

A - Po^o/Z (34) 

where Z/pqCq is the nondimens ional normal incidence impedance. 

The solution for acoustic pressure within the cavity is found from a 
weighted residual statement: Find a function p(x,y,z) from the class of con- 

tinuous functions which satisfies (Vp are the values of Vp on the boundary) 

; W[v2p. (1/cJ) Ptt)d« - / W[vp-vp].n dS = 0 (35) 

Q S 

/s 

for all test functions W from the class of continuous functions, where 3Q 

takes on the appropriate value on the constituent segments 
of 30, and Vp are the values of Vp on the boundary, from Eq. (32). 

By using the Divergence Theorem twice on the volume integral, and by 
defining members of the set of test functions to be solutions of the eigen- 
value problem 

v2^ + (w/c)2^ - 0 (36) 

V^*n = 0 on dO 

the weighted residuals statement becomes 

f ^itPtt + w^^p]dn = -Po^o^ f -Cq X Pt dS (37) 

“ L 
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where is an eigensolution to Eq.(36) corresponding to the eigenvalue co^. 
These eigensolutions are defined by Eq.(30). 

The function p(x,y,z) is expanded in terms of the complete sequence of 
functions so that 

p(x,y,z) « [^]{a) 

where [<f>] is a row matrix whose elements are a truncated sequence of the func 
tions defined by Eq. (36) and (a) is a column matrix of time -dependent 
amplitude coefficients. 

The acceleration on the sidewall is similarly expanded in terms of the 
structural eigenfunctions evaluated on the plate surface. 

Wtt(x.z) - [Np] [^p] (q) 

where the row matrix [Np] is the implicitly- defined shape matrix for the side 
wall deflections and [V>p] is the sidewall partition of the modal matrix 
derived from the structural free vibration problem. The vector (q) is a vec- 
tor of time dependent amplitude coefficients for the structural modes. Equa- 
tion (32) becomes 


where 



-[D]{q) 


(38) 



(39) 


'bO. 

s 

The diagonal matrix is the "generalized mass matrix" for the cavity 

(actual units are volume), [C] is generally not diagonal and is the general- 
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ized damping matrix (the form here is only valid for harmonic motion since the 
impedance relation [28] only applies for harmonic motion). The matrix [D] is 
the coupling matrix and is generally not square. If N^, cavity modes and Ng 
structural modes are retained in the model, then [D] is Nj, x Ng. Note that 
the matrix [D] which couples the plate into the cavity, appears as [D]^ in the 
coupling of the cavity into the plate. 


THE PROPELLER GENERALIZED FORCE 

The conclusion to be drawn from Eq. (17) is that the lift distribution 
(lift per unit span) induced at a point on the wing can be written in the form 

^ 0 ) ~ f ^(>?)e^^N^^ 

where ^ 

L - (l/2)pu2cC (r/27TRU) 

A/ 

a 

w 

and the nondimens ional spanwise distribution functions f^(y) are defined over 
four regions by Eqs . (1) through (4). 

The generalized force is determined by the virtual work integral as in the 
determination of the acoustic pressure generalized force in Eq.(28). In the 
case of the propeller loading 

s 

=* ^ 5w(x)dx 
o 

where x is the wing spanwise coordinate and S is the span. The wing vertical 
deflection is given in terms of the implicitly defined global shape matrix 
[N^] as 

w(x) = [N^]{W^} 

where (W^) is here defined to include the nodal displacements and bending 
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slope at the attachment point. Then 


5W - 


£3 

^ [N^]^ 


(x)dx 


Consistent with finite element procedures, the integral is carried out over 
the element subdomains and assembled into the global generalized force so that 


5W 


■^^{5W^)ny' [N„®]T je^(x)dx 


n 


n 


where the integral is over the n^^ subdomain. The element shape function 
matrix is [N^®] . The element generalized force vector is thus 


w 


where 


il n 

-^(QwW 

s, 


tQw^ni “ ^ y f^(r;)dx 

n 

It follows that the global generalized force vector is of the form 


(40) 


(41) 


(42) 


-s 


i 


iiNQt 


(43) 


where (Q„) is obtained by assembling the {Qw)ni conventional methods. 


THE AIRBORNE GENERALIZED FORCE 

The pressure level in the acoustic field of the propeller is given by 


p(x 


,r,0,t) ^ P^(x,r)e^ 


(S,NS)t-£N0) 


£=--■ 
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The finite element analysis yields the pressure field in each angular har- 
monic. That is, we obtain Pjg(x,r) corresponding to the angular harmonic iN 
which is related to the J0th harmonic of blade passage frequency iNO. 

As shown in Figure 4, the acoustic field of the propeller at a fixed 
radial distance from the propeller axis shows considerable variation with 
axial location, that is the propeller radiated field varies for and aft. In 
the example shown in the figure, the maximum sound pressure level occurs aft 
of the propeller and is about 13 dB above a second relative maximum which 
occurs ahead of the propeller. 

At the present time little is known about the coupling of the propeller 
acoustic field and an adjacent structure. The field is neither normally inci- 
dent nor representable by oblique plane waves in a simple way. In this study 
we have used the simplest estimate of the driving pressure due to the propel- 
ler noise field. The level is taken as uniform on the sidewall at some frac- 
tion of the peak level. The results actually reported here use the peak level 
itself. At lower fractions of the peak level the interior noise results can 
be expected to be reduced accordingly. 

With this assumption the sidewall is subjected to the distributed normal 
load 

f(x,z) “ ^ giiNQt 



where is the pressure amplitude in the harmonic of the blade passage 
frequency chosen, as noted above. The virtual work expression becomes 


5W=-{fiWp)y* [Np]T^] 

(L 


Pjl 


At the element level the generalized force is 


(44) 
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/ _ , iJlNfit 

[Np]T (P_j) dS e (45) 


and after assembly the system generalized force is of the form 


(Qs) °y^(Qs)^ e^^NQt (46) 

I 

which is again a superposition of harmonics of the blade passage frequency. 


THE COUPLED EQUATIONS OF MOTION 

The system equations of motion can now be written. For the structural 
system 





• 


(q) + 

i^si^sii 

(q) + 

“si ^sii 


- ^ 




(q) 


[D]T{a) » (Qi) 


(47) 


For the acoustic cavity 


“1 

,, • 


^cii 
_ ^ 

(a) + [c]{a) + 

_ X 


(a) + = (0) 


(48) 


The structural system consists of Ng equations and the cavity system consists 
of equations. Equations (47) and (48) thus represented a system of (Ng + 
N^.) equations set up for dynamic response calculated with the input vector 

(Qi) - [^]^{Qw> 

or (49) 

(Ql) = [\&]'^{Qs) 

Reference to equations (43) and (46) show that {Qj) can be written in terms of 
its frequency components 

(Ql) -z ei^NOt (50) 

a 
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where for structureborne inputs 

{Qi)i - [^]'^{Qw)i 

and for airborne inputs 

COMPUTATIONAL RESULTS 

In this report, results are presented for the comparison of interior lev- 
els for structureborne and airborne noise for several cases of the structural 
stiffness and mass parameters. For this purpose we have chosen a specific 
propeller model with the following characteristics: 

I 

Number of blades - 4 

Radius - 7 ft. (2.133 m) 

Blade chord (constant) - 1 ft. (.305 m) 

Blade section lift coefficient - 2n 

Thrust - 4000 lb (17800 N) 

RPM Range - 750-1500 RPM 

Tip-Fuselage clearance - 7 ft. (2.133 m) 

Flight Mach number - 0.27 

Speed of sound - 1125 ft/sec. (343 m/sec.) 

The thrust is held constant for the range of propeller speeds so that the 
blade twist and loading varies with RPM for this study. This also means that 
the top loading and therefore the trailing vortex strength varies with RPM. 

The wing model employed is one of span 29.5 ft (9 m) and chord 7 ft. 
(2.133 m) with a section lift coefficient of 2n, The propeller axis is 
located 18 ft. (5.5 m) from the wing root. The reference point for specifi- 
cation of the airborne noise level is at 14 ft. (4.27 m) from the propeller 
axis. This is a somewhat arbitrary choice in that it places the effective 
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wing root inside the fuselage sidewall, however, only minor changes in inte- 
rior levels due to the airborne noise are associated with this choice. 

Figure 7 shows the computed maximum sound pressure levels in the radiated 
field of the propeller at two propeller radii from the axis. Levels are shown 
for the first three harmonics over a propeller speed range of 750 - 1500 RPM. 
For the propeller chosen, the maximum speed is close to sonic tip velocity. 
Particularly at the lower rotational speeds the sound pressure level drops off 
rapidly with increasing harmonic number. The levels shown in Figure 7 are 
used as the input for the computation of the airborne contribution to interior 
noise . 

The interior cavity is 19.7 ft (6 m) in length, 6.88 ft (2.1m) in height, 
and 7.87 ft (2.4 ra) in depth. The acoustic characteristics of the cavity are 
such that 80 acoustic modes cover the frequency range up to 320 Hz. 

The structural properties of importance are the mass and stiffness of the 
wing and carry through structure, and the mass and stiffness properties of the 
sidewall. All results presented here are for one sidewall configuration. The 
sidewall is taken as an aluminum plate of thickness approximately .080 in. 
(.002 m) . Relatively light frames and stringers are treated as smeared mass 
and stiffness yielding a stiffness properties matrix 


where 



(51) 


= D + (Eglg/d) 


= 


= Eh^/[12(l-u2) ] + (Egls/d) “ 16.4x10^ Ib-ft (22.3x10^ N-m) 
D + (Eflf/h = 10.2x10^ ft-lb (13.9x10^ N-m) 

Dvv = uD = 12.9 ft-lb (17.5 N-m) 
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[(l-v)/2] D - 13.1 ft-lb (17.8 N-m) 


The smeared mass of the plate is 

^e “ ^ (PfAf/^) - 0.054 slugs/ft^ (8.55 kg/m^) 

In Eqs . (51) and (52), Eg and Ef are the Young's modulus for the stiffeners 

and frames, Ig and If are the cross-sectional area moments of inertia of the 
stiffeners and frames, d and i are the element edge lengths corresponding to 
the stiffeners and frames, pg and pf are the mass per unit length of the 

stiffeners and frames, Ag and Af are the cross-sectional areas of the stiffen- 
ers and frames, and Mp is the plate mass per unit area. The plate Young's 

modulus is E and the plate thickness is h. Eight cases of wing and carry- 
through structure stiffness and mass have been considered: 

Case I : Stiff Light Wing and Light Carry Through 

Wing: El = 10^ Ib-ft^ (0.414 x 10^ N-m^) 

p « 1,13 slugs/ft (54 kg/ra) 

Carry Through: El = 10^ Ib-ft^ (0.414 x 10^ N-m^) 

p ^ 1.13 slugs/ft (54 kg/m) 

Wing Fundamental Frequency — 19 Hz 


Case II : Flexible, Heavy Wing and Stiff, Light Carry Through 

Wing: El = 0.5 X 10^ Ib-ft^ (0.207 x 10^ N-m^) 

p =* 5.63 slugs/ft (270 kg/m) 

Carry Through: El = 10^ Ib-ft^ (0.414 x 10^ N-m^) 

p = 1.13 slugs/ft (54 kg/m) 

Wing Fundamental Frequency = 6 Hz . 
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Case III : Flexible, Heavy Wing and Flexible, Heavy Carry Through 

Wing: El - 0,5 X 10^ Ib-ft^ (0,207 x 10^ N-m^) 

p - 5,63 slugs/ft (270 kg/m) 

Carry Through: El - 0,5 x 10^ Ib-ft^ (0.207 x 10^ N-m^) 

p - 5.63 slugs/ft (270) kg/m) 

Wing Fundamental Frequency - 0.2 Hz. 

Case IV : Flexible, Heavy Wing and Flexible, Light Carry Through 

Wing: El - 0.5 X 10^ Ib-ft^ (0.207 x 10^ N-m^) 

p = 5.63 slugs/ft (270 kg/m) 

Carry Through: El - 10^ Ib-ft^ (0.414 x 10® N-m^) 

/9 = 1.13 slugs/ft (54 kg/m) 

Wing Fundamental Frequency - 0.2 Hz. 

Case V : Case I with Acoustic Damping 

A = 0.5, where A is defined by equation (34). 

Case VI : Realistic Wing and Realistic Carry Through 

Wing: El = 10® Ib-ft^ (0.413 x 10® N-m^) 

p = 4.532 slugs/ft (217 kg/m) 

Carry Through: El = 10® Ib-ft^ (0.413 x 10® N-m^) 

p = 4.532 slugs/ft (217 kg/m) 

Wing Fundamental Frequency = 3 Hz . 

Case VII : Realistic Wing and More Flexible Carry Through 

Wing: El = 10® Ib-ft^ (0.413 x 10® N-m^) 

p = 4.532 slugs/ft (217 kg/m) 
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Carry Through: El - 0.5 x 10® lb- ft (0.206 x 10® N-m) 

p - 4.532 slugs/ft (217 kg/m) 

Wing Fundamental Frequency — 3 Hz . 

Case VIII : Realistic Wing and Stiffer Carry Through 

Wing: El - 10® Ib-ft^ (0.413 x 10® N-m^) 

p = 4.532 slugs/ft (217 kg/m) 

Carry Through: El = 2 x 10® Ib-ft (0.826 x 10® N-m) 

p * 4.532 slugs/ft (217 kg/m) 

Wing Fundamental Frequency ^ 3 Hz . 

Case I represents a very stiff and relatively light wing and carry- through 
structure. Case II is a somewhat less stiff, but significantly more massive 
wing structure with the carry- through structure of Case I. Cases III and IV 
are at the other extreme. Both have a flexible, relatively massive wing 
structure. In Case III the carry- through structure is flexible and massive 
while in Case IV it is flexible and relatively light. While not conforming to 
any existing structures, these cases should at least bound realistic config- 
urations. Case V is just Case I with interior acoustic damping. Cases VI, 

VII and VIII are much closer to realistic structures. Case VI establishes a 
representative wing and carry- through system. In Case VII, the stiffness of 
the carry- through structure is one-half of that of Case VI; in Case VIII, it 
is twice as stiff as in Case VI. 

The structural model for the dynamic response calculations is based on a 
subset of the structural eigenvalues and eigenvectors spanning a frequency 
range in excess of 1000 Hz, In the present study we have not included struc- 
tural damping and we have not included dissipation within the cavity. Struc- 
tural damping has the primary effect of limiting the resonant response. It is 
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not quantifiable with any degree of certainty. 

Any effort to assess the interior noise levels on the basis of resonant 
response amplitude is not appropriate due to the sensitivity of the response 
amplitude to damping. Our assessment of the relative levels of interior noise 
from the airborne and structureborne sources is based on a comparison of the 
non-resonant response, that is, a comparison of the general levels of the 
response . 

For both structureborne and airborne inputs, equation (50) shows that the 
generalized force is a superposition of harmonics of the fundamental frequency 
NO. In the present study we retain three harmonics since it has been found 
that the nonresonant response levels are reasonably well defined with this 
limited representation of the input. The interior response will also exhibit 
the harmonics of the fundamental. The overall acoustic pressure level can 
therefore be represented simply in terms of the amplitudes of the harmonics 
according to 

T r 3 

p2 - ly* I p|2 dt - I 

^ i=l 

where p^ is the amplitude of the pressure response at the measurement point 
for the i^^ harmonic of the fundamental frequency NQ and p^ is the average 
square amplitude of the response pressure. 

The results of calculations for Cases I through VIII are shown in Figures 
8-25. These figures show the interior noise level at a point adjacent to the 
fuselage side wall in the center of the sidewall. The response is shown as a 
function of propeller rotational speed in RPM. At each RPM the input has many 
harmonics (three have been considered here) so that the noise level is defined 
by equation (53) superimposing the effects of the harmonics. On Figures 8-11, 
the structureborne and airborne levels are shown separately. Note that the 
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two curves have some resonant peaks in common and others at different 
frequencies. This is because of the different nature of the structureborne 
and airborne generalized forces and the modes (both cavity and structural) 
which they excite. Shown superimposed on the response curves are interpreta- 
tions of what we define as the general level of the nonresonant response which 
because of the uncertainty in structural damping level is used to assess the 
comparison between the two noise sources. 

Figure 8 shows the results for Case I in which the wing and carry- through 
structure is very stiff. The wing cantilevered fundamental frequency is 19 
Hz. It is noted in this case that the structureborne noise creates levels 
substantially below those related to airborne noise. 

In Case II, shown on Figure 9, the wing is considerably less stiff than in 
Case I. The cantilevered fundamental frequency is 6 Hz. The carry- through 
structure is still extremely stiff and light. The interior levels in this 
case are still dominated by the airborne noise. 

In the third and fourth cases, the wing and carry- through structure are 
very flexible. The cantilevered wing fundamental frequency is 0.2 Hz in both 
cases, perhaps a factor of ten below actual wing characteristics. In Case III 
the carry- through structure is less stiff than in Case IV. Figures 10 and 11 
show that the interior noise level contributions of the structureborne and 
airborne sources become much more similar, particularly at the lower propeller 
speeds where the two levels are about the same. 

It is clear that the stiffness and mass of the wing and carry- through 
structure are of central importance in the transmission of structureborne 
noise into the fuselage. In all four cases the airborne noise level is rela- 
tively unaffected by the wing and carry- through structural characteristics. 
However, the structureborne noise level is strongly dependent on these char- 
acteristics. The structureborne levels generally decrease with an increase in 
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the stiffness of the wing and carry-through structure. 

Figures 12 and 13 show the effect of interior acoustic damping on airborne 
and structureborne noise. In both cases, damping (as expected) is seen to 
damp out peaks. Especially affected are those peaks that are due to acoustic 
resonances . 

Figures 14-16 are for the more realistic wing and carry- through structure. 
For these cases, airborne noise is significantly higher than structureborne 
noise. These results are consistent with the previous cases in the sense that 
lower structureborne noise levels are associated with increased stiffnesses of 
the carry- through structure. 

Figures 17-19 show airborne noise levels when the fuselage is only par- 
tially loaded. In these calculations, the cases VI, VII, and VIII are 
repeated with only the panel loaded in a centrally- located vertical strip 
covering 60% of panel area. Both damped and undamped cases have been com- 
puted. Comparison with Figures 14-16 indicates that there is only an insignif- 
icant change in the airborne noise levels when the panel is only partially 
loaded. 

In Figures 20-25, the effect of various damping levels is shown on the 
airborne and structureborne noise levels for cases VI-VIII. As expected, the 
damping reduces peaks in the response that result from cavity and structure 
resonances. The meaning of A is given by equation (34). 


CONCLUSIONS 

A simple, though physically meaningful, model for assessing the relative 
importance of structureborne and airborne contributions in interior noise 
levels for propeller driven aircraft has been created. The system for which 
dynamic response calculations are made consists of a cavity model, a struc- 
tural model, a model for the propeller radiated acoustic field, and a model 
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for the propeller trailing vortex mechanism which excites the structureborne 
noise. The dynamic response model is based on the cavity modes and frequen- 
cies and the wing - carry through structure - sidewall modes and frequencies 
with suitable coupling. The forcing function is synthesized from the gener- 
alized forces associated with the propeller radiated field and the trailing 
vortex system. 

Computational results indicate that the relative importance of the two 
source mechanisms depend strongly on the stiffness and mass properties of the 
wing carry- through structure. For stiff structures the structureborne inte- 
rior noise levels are well below the airborne levels. For very flexible 
structures the interior levels from the two source mechanisms can be compara- 
ble. 

Both the geometry of the structural model and the stiffness and mass 
properties have been changed for the purpose of investigating different con- 
figurations with the goal of more effectively bracketing the characteristics 
encountered in actual aircraft installations. The structural model can employ 
structural damping or equivalent viscous damping. At this stage we have cho- 
sen not to use it because of lack of quantifiable data. The cavity model did 
include dissipation due to acoustic treatment. The effect of internal 
absorption was studied. 

The noise source models are constructed so that simple data changes allow 
for the variation of all important propeller parameters. It is a simple 
matter to investigate propellers of any radius and any number of blades, 
positioned as required relative to the fuselage. 
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A typical propeller radiation pattern. Contours of constant acoustic 
pressure magnitude. The propeller has six blades and the tip speed is 
approximately sonic. 






Figure 5. The idealized structural model with the fuselage cavity 
modelled by a rectangular cavity with a flexible side 
wall. The wing is modelled as a beam carried into the 
fuselage by a discrete sitffener. 
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Figure 6. The finite element model of the fuselage side wall 
showing the element subdivisions and the discrete 
stiffener for the carry- through structure. The 
elements can have smeared or discrete stiffeners. 
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Figure 8. Interior noise levels for structural model. Case I, as a function 
of propeller RPM at constant thrust. 
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Figure 9. Interior noise levels for structural model. Case II, as 
function of propeller RPM at constant thrust. 




200 



Q 

in 



COP) |;aAa“i asio/g jorja-^uj 


43 



STRUCTURAL 
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Interior noise levels for structural model. Case IV, as 
function of propeller RPM at constant thrust. 



200 


I 

I 



45 


Figure 12, Effect of intrior absorption on interior noise levels due to 
airborne noise transmission for the structural model of Case 
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structurborne noise transmission for the structural model of Case 
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STRUCTURAL AND AIRBORNE NOISE TRANSMISSION (CASE VII) 
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STRUCTURAL AND AIRBORNE NOISE TRANSMISSION (CASE VIII) 
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Figure 17. Interior noise level due to airborne transmission, partial 
loadi ng for Case VI . 
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Figure 18. Interior noise level due to airborne transmission, partial 
loadi ng for Case VII . 
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